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We review the self-consistent mean-field theory for charge-frustrated Josephson junction arrays. 
Using (cos (fi) {(f is the phase of the superconducting wavefunction) as order parameter and imposing 
the self-consistency condition, we compute the phase boundary line between the superconducting 
region ((cos<^) 7^ 0) and the insulating one {{cos(f) — 0). For a uniform offset charge q = e the 
superconducting phase increases with respect to the situation in which q = 0. Here, we generalize 
the self-consistent mean-field theory to include the effects induced by a random distribution of 
offset charges and/or of diagonal self-capacitances. For most of the phase diagram, our results agree 
with the outcomes of Quantum Monte Carlo simulations as well as with previous studies using the 
path-integral approach. 



I. INTRODUCTION 

The first artificially fabricated Josephson junction arrays (JJA's) were realized twenty years ago at IBM Q as an 
effort to develop an electronics based on superconducting devices. Immediately after, it became clear that JJA's 
provided an ideal model to investigate classical phase transitions, frustration effects and relevant aspects of non-linear 
dynamics 0, JJA's are built by placing, on the sites of a lattice, islands of superconducting material coupled 
by Josephson junctions. The huge variety of behaviors of the system is rather simply described by the competition 
between the Josephson energy Ej and the charging energy Ec' the former being responsible for the Josephson 
tunneling of Cooper pairs between the sites of the lattice, while the latter measures the effects of the electrostatic 
repulsion between Cooper pairs. The superconductor-insulator transition typical of JJA's, for instance, depends 
crucially on the ratio between these two energy scales. 

In many situations, it is relevant to analyze the effect of a background of external charges on the superconductor- 
insulator transition of a quantum JJA. Offset charges arise in experimental realizations of JJA's as a result of charged 
impurities or by the application of a gate voltage between the array and the ground. In the former situation, offset 
charges are naturally randomly distributed on the lattice while in the latter situation they play the role of a sort of 
chemical potential and, then, their distribution may also be uniform. Offset charges may be regarded as effective 
charges qi located at the sites of the lattice: when qi ^ 2e, they cannot be eliminated by Cooper pair tunneling. 
In general, offset charges frustrate the attempts of the system to minimize the energy of the charge distribution of 
the ground state (for this reason they are also called frustration charges). A large number of studies has by now 
been devoted to the analysis of the effects induced by offset char ges both on the zero-temperature phase transition 
^lEIEQ and on the phase transition at finite temperature |^l^[lof. 

In this paper we shall use the self-consistent mean-field theory (SCMFT) to investigate the finite temperature phase 
diagram for the self-charging (SC) model of JJA's 0; in this model it is assumed that the potential at site i depends 
only on the charge at the same site and, thus, the capacitance matrix describing the charge effects of the array is 
diagonal. In particular, we shall investigate situations in which offset charges (both uniform and random) are present. 
Although quantum corrections may be relevant for d < 2, the SCMFT has the merit of providing a rather intuitive and 
physically transparent approach to the analysis of some general features of the superconductor-insulator transition in 
these systems. Furthermore, the results we obtained are in very good agreement with the outcome of recent numerical 
simulations and consistent with other analytical approaches not relying on mean- field theory [T^ . 

The plan of the paper is as follows: in Sec. II we review the SCMFT for the SC model of quantum JJA's with a 
uniform distribution of offset charges Il0j|. We study the eigenvalue equation of the mean-field Hamiltonian and, 
for a uniform offset charge qi — e, we show that there is superconductivity for all values of Ej/Ec. In Sec. Ill we 
point out how to extend the SCMFT to situations in which capacitive disorder is present: one has to impose the 
self-consistency condition with a double average, the quantum one and the average over the disorder. The results 
are in agreement with the ones obtained with the path-integral approach and, at very low temperatures, are 
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consistent with the phase diagram obtained in Ref. [l5|- Section IV is devoted to our concluding remarks. 

II. MEAN-FIELD THEORY FOR JJA'S WITH OFFSET CHARGES 

The Hamiltonian commonly used to describe Cooper pairs tunneling in superconducting quantum networks defines 
the so-called quantum phase model (QPM): 

y (ij) 

where ipi is the phase of the superconducting order parameter at the grain i. Its conjugate variable rii = Qi/2e — 
—id/dLpi (with [(/3i,rtj] = i Sij) describes the number of excess Cooper pairs on the i-th superconducting grain and 
Cij is the capacitance matrix. The symbol (ij) denotes a sum over nearest-neighbor grains only. 

The first term in the Hamiltonian ^ determines the electrostatic coupling between the Cooper pairs while the 
second term describes the hopping of Cooper pairs between neighboring sites (Ej is the Josephson energy). An 
external gate voltage Vi provides a contribution to the energy via the offset charge qi = Cij Vj ; this external 
voltage can be either applied to the ground plane or it may be induced by charges trapped in the substrate. The 
former situation leads to the appearance of a uniform frustration charge, while the latter naturally induces a random 
offset charge. In this paper we shall limit our investigation only to the SC model described by the Hamiltonian 
with Cij = 6ij Cii. When all the Ca^s are equal {Ca — Cq), the charging energy Ec is defined as Ec = /2Cq. 

With a uniform distribution of offset charges qi = q, the Hamiltonian of the array is given by: 

H = AEcY.{^^^-^^^ -E.jY,cos{^,^^,). (2) 



Mean-field theory for the SC model for quantum JJA's was first used by Simanek 0, IT^ . The approximation 
consists in replacing the Josephson coupling on the link i-j by an average coupling so that Ej J2{ij) cos{ipi — ipj) « 
zEj{cos (f) cos ifi, where z is the coordination number. Requiring the order parameter to be real, leads to (sin (ft) = 
0; it is also assumed that (cos 93) does not depend on the island index i. In the mean-field approximation the 
Hamiltonian (j^Jis given by a sum of single site Hamiltonians Hi describing a quantum particle in the potential cos ipi: 

The pertinent Schrodinger equation to be solved is then 

-AEc-r^ - &i-^Ec-r- + 4 (■^) Ec ~ zJ{cosip) cos J V'm((^) = E,nTpmi(p)- (4) 
L clcp^ Ze dip \IeJ J 

Due to the periodicity of the phase the eigenfunctions should be 27r-periodic functions of (f, i.e., 

4>n {^) = i'n (<P + 27r) . (5) 

Furthermore, since the Hamiltonian is invariant under the shift q ^ q~\- 2Ze, where / is an integer, it is relevant to 
analyze only the situations corresponding to g = and q ~ e. 

The order parameter (coscp) is evaluated in terms of the eigenfunctions of Eq. through the self-consistency 
equation 

Ee"''"/'"^(V'm|cos(/3|V'™) 
= ^e-^--./fc^^ • 

ni 

From Eq. JBJ, one immediately sees that, for high temperatures or low Ej, only the solution {cos ip) = exists 
and, thus, there is no superconductivity; for low temperatures or high Ej instead, (cos(/?) and the system as 
a whole behaves as a superconductor. Solving the eigenvalue equation 10} provides us with all the tools needed to 
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investigate the finite temperature phase diagram of the SC model of frustrated JJA's. Defining v = —zEj{cos (p) /2Ec', 
A™ = [E^ - 4:Ec{q/2e)'^]/Ec, and K = 2i{q/2e), one finds 

Equation ((TI) is a Mathieu equation with a term proportional to a first derivative: setting 

^™((^) =e-^^VmM, (8) 

one gets an equation for pm, namely 



If one sets A„i = Ep^/Ec and tp = 2x, the eigenvalue equation jT)) is usefully recast in the standard form of the 
Mathieu equation |17|: 

+ (A„ - 2v cos 2x)pm = 0. (10) 

It is well known that the Mathieu equation admits the following periodic solutions: 
1- ce2n{x,v), even solutions with period tt corresponding to the eigenvalues a2n{v)', 

2. se2n+2{x,v), odd solutions with period tt corresponding to the eigenvalues b2n+2{i^)] 

3. ce27i+i{x , v) , even solutions with period 27r corresponding to the eigenvalues a2n+i{v)] 

4. se2n+iix,v), odd solutions with period 27r corresponding to the eigenvalues 62n+i(*^)- 

Since the eigenf unctions of the Schrodinger equation |0J should satisfy the periodic boundary condition jS)), from Eq. 
^ one immediately sees that one should treat differently the situations where q/2e is integer or half-integer. In fact, 
for integer q/2e, one has to consider only 7r-periodic solution of Eq. (|10(l : in this way Pmip) are 27r-periodic, which, 
in turn, leads to 27r-periodic ipmif) — Pmip) (since / is an integer). The solutions of Eq. 1)10(1 with period tt are 
the Mathieu eigenfunctions ce2n+2{x), se2„(x) (with n = 0, 1, . . .) |l7l |. If, instead, q/2e is half-integer, Eqs. jSJ) and 
© require the use of the Mathieu eigenfunctions ce2„+i(a;), se2rt+i(a^) (with n = 0, 1, . . .) which are tt- anti-periodic. 
Then Pmip) are 27r-anti-periodic and the eigenfunctions V'(^) = Pm{p) are 27r-periodic. 

Since the phase transition is expected to be second order |lflj . near the transition temperature, the order parameter 
{cos if) and the parameter v are small: this allows one to use the expansion for small v of the Mathieu functions [l7l |. 
As a result, one finds that, at the first order in v and apart from the phase factor e^^'^^^, the normalized eigenfunctions 
of the Schrodinger equation satisfying, for q — e, the condition Q are given by 



If V 3(p\ 

cos 

2 8 2 J 



V;i(^) = -=(sm---sm- 



(11) 



and, for n = 1, 2, . . ., by 

V'2„+i(^) = 2 ^ 



/r. , i\ (2n+3)ip (2n-l)<fi 

{2n+l)Lp ^cos^ — cos^^ — 



4(2n + 2) 8n 



1 /r^ , 1 \ ■ (2n+3)uP • (2n—l)ip 

1 r. {2n+l)ip [sm ^ ^ >^ sm^ ^ 



(12) 



V 



4(2n + 2) 8n 

where the superscript e (o) stands for even (odd). The corresponding energy eigenvalues are 

E1^Ecil + v)/2, 
El^Ec{l-v)/2, 

£^2„+i=^2Vi-^c(2n+ 1)2/2 
n = l,2,... 



(13) 
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The expectation values (V'ml cos Lp\'4)m), at the order v, are given by 

1 V 

cos if = 2^8 



(V'2n+l|cOSV?|?/'2„+i) = (V'2n+l|cOS</3|V'2„+l) 



8n(n + 1) 
n = l,2,... . 

Upon inserting the above eigenfunctions and eigenvalues in Eq. (jBJ and keeping only the terms proportional to 
V ^ (cosiys), one obtains the following equation for the critical temperature Tj,: 

-=5(g = e,y), (15) 
a 



where y = ksTc/Ec, a = zJ/{AEc) and 



Ay ^ l-4(n+l/2)2' 

5(9 = e,y) = "^^ 00 • (16) 



-1/y + ^ g-(4/y)(«+i/2)2 



=1 

From Eq. H15(l one easily shows that, in the presence of charge frustration ±(2n + 1) e on the lattice sites, for each 
value of a, there is a insulator-superconductor transition. Indeed, g{e, y) — > 00 for y and g{e, y) — > for y — > 00; 
also dg/dy < for all y > 0. It follows that Eq. 115|) has a unique solution for each value of a. Moreover, since 
g ~ 1/y — 1/a, the critical temperature at which the transition occurs is given by Tc ~ zEj/A. 

For q — 0, the solutions of the Mathieu equation (|10|l with period tt are the Mathieu eigenfunctions ce2n+2{x), 
se2n{x) (with n = 0, 1, . . .): one finds 

--.g(g = 0,y) (17) 
a 

where 

5(9 = 0,y) = . (18) 

1 + 2 e~4»Vy 

n=l 

From Eq. p7|l one sees that for a > 1 the system has an insulator-superconductor transition, while for a < 1 there is 
no evidence for a transition. In fact, g{e, y) — *■ 00 for y — > and g{e, y) — > 1 for y — ^ cxd; also dg/dy < for all y > 0. 
Therefore, the self-consistency equation © does not have solutions for a < 1 and it has a unique solution for a > 1. 
For small values of y, Eq. ^ yields 1/a « (1 - 2e-^/y/3)/{l + 2e-^^y), from which 

ksTc 4 

^ i^^T^Is+lI' ^^^^ 

^"6 3(a-l) 

In Fig. ^we plot Tc as a function of a for (7 = and g = e using Eqs. p5|) and (|17|) . For q — 0, one sees that there 
is no superconductivity for a < 1. For q — e, superconductivity is attained for all the values of a: a uniform offset 
charge q = e always favors superconductivity. 

It is worth noting that the eigenfunctions ipln' of the Schrodinger equation without the periodic potential are 

^l^'> ^ ^e^™'^/^. (20) 

These wavefunctions are also eigenfunctions of the number operator with eigenvalues J\f = {ipm^l — id/dip\ipin) = 
±m/2. When q — 0, the wave functions V-'in and "0271 are, respectively, the even and odd combinations resulting 
from the splitting of the eigenfunctions ip^J and are related to the expectation value of the half-integer number of 
Cooper pairs {J\f = 0,±1,±2, . . .). On the other hand, when q = e, the expectation value of the number operator 
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FIG. 1: Phase diagram for the SC model without charge frustration (sohd hne) and with charge frustration g = e (dashed hne). 
a stands for the ratio zEj/AEc- The I and S indicate, respectively, insulating and superconducting phase. 



on the eigenfunctions 'il'2n+i ^^'^ V'fn+i half-integer and is equal to 1/2 in the ground state. Therefore, an offset 
charge q ~ e favors the Cooper pairs tunneling, making possible the insulator-superconductor transition also when 
Ec > Ej. 

If one should use both periodic and anti-periodic solutions, the general solution of Eq. (@J) would not have a 
definite periodicity and, consequently, the charges would take any value; this situation is expected to be relevant 
in the description of continuous flows of current due, for instance, to ohmic shunt resistances 113,01. Unless there 
is dissipation, the use of both periodic and anti-periodic solution is unwarranted; however if, in the self-consistency 
equation ©, one should include also the 27r-anti-periodic eigenfunctions 'IG'I, one would find - for small critical 
temperatures - the following equation 

1 = « ^T7 T7 • (21) 

1 + 2e-^/y + 2e-4/y ^ ' 

Equation (|21|1 for a less than a critical value ac (etc ~ 0.79) does not have solutions, for ac < a < 1 lias two solutions 
and for a > 1 has just one solution. This behavior is called reentrant [l6|. 

We conclude this Section observing that the phase boundary line obtained within the mean-field approximation in 
the path-integral approach is 1/a — g{q, y) with J), 10] 

5(9, y) = ^ -i(^+jL.. (22) 



e 



Equation 1221) coincides with Eqs. H15|) and (|17|l for g e and (7 = 0, respectively. 



III. CAPACITIVE DISORDER 

In practical realizations of Josephson devices 0, one has to deal with disorder caused by offset charge defects 
in the junctions or in the substrate |23|. Random offset charges cannot be made to vanish by using a gate for 
each superconducting island since, in large arrays, too many electrodes would be necessary, making impossible the 
cooling of the system at the desired temperatures. In Ref. [2J| it was observed a sensible variation (~ 40 %) of the 
resistance between the unfrustrated and the fully frustrated array. Moreover, it may also happen that the network's 
parameters are not uniform across the whole array: despite recent advances in fabrication techniques, variation of 
junction parameters associated to the shape of the islands can be also of 20% % Thus, it is relevant in many practical 
situations to study JJA's also with randomly distributed self-capacitances; this corresponds to have a random diagonal 
charging energy 0, . 

In this Section, we shall determine the finite temperature phase diagram of JJA's with capacitive disorder (i.e., with 
random offset charges and/or random self-capacitances). To derive the phase boundary between the insulating and 
the superconducting phase, we shall use the mean-field approach for quantum JJA's with offset charges and diagonal 
capacitance matrices reviewed in the previous Section. One has to impose now the self-consistency condition with a 
double average: the quantum average and the one over the disorder. 

As we shall see, charge disorder supports superconductivity; furthermore, the relative changes of the insulating 
and superconducting regions of the phase diagram depend crucially on the weights of the (5-like charge probability 
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distribution. In the physical relevant situation of two charge distributions peaked at the values q = and q = e, 
increasing the frustrated weight favors the superconducting phase. Also the randomness of the self-capacitances leads 
to remarkable effects, namely, the superconducting phase increases with respect to the case where disorder is not 
present. In the following, we shall provide a quantitative analysis of these phenomena. 

A pertinent extension of SCMFT in the presence of on-site disorder may be obtained introducing an order parameter 
averaged also over the disorder. In the following (• • • ) denotes only the quantum average while [• • • ]av an average over 
the random variables. The single-site Hamiltonians of Eq. (|3Jl then become 

H^ = -'^^Cg^ - ^'^^'^£9^ " [(C0S^)],„C0S(P. (23) 

where i^^-* — e'^C'^^ /2. The Hamiltonian (|23|l depends on a random variable X, which can be cither qi or E^} . Thus, 
its eigenfunctions and eigenvalues depend either on qi or E^^ : 

H,^n{w,X) = En{X)M^,-X). (24) 

The self-consistency condition is given by 

[/ M A.^ o.^^ E„e-^^"'^^(V'n(^)|cos^|V^„(X)) 

[(cos = J (X) j2 e-PEAX) (25) 

where P (X) is the probability distribution of X. 

The phase boundary line is obtained from Eq. (|25|l by requiring [{cos ip)]^^ to be small and by keeping only terms 
proportional to it. The self-consistency condition yields a mean- field phase boundary line in agreement with the results 
obtained by the path- integral approach [T^ . The low temperature behavior, obtained by a pertinent extrapolation of 



rjObtamed by ape 

our finite T results, is consistent with the phase diagram obtained in Ref. [l5| (see Ref. [ij) 



A. Random Offset Charges 

In this Section we shall study JJA's at finite temperature with random charge frustration. One assumes that 
the offset charges are independently distributed according to a probability distribution given by a sum of (5-like 
distributions 

n 

with X^nP'i ~ 1- This corresponds to a random distribution of charges which are integer multiples of e and, actually, 
this is the most realistic situation for a random distribution. Inserting the probability distribution (|26() in Eq. I|25|) . 
one has 

~ = J dq^PnSiq-ne)g{q,y) ^^png{ne,y)+^png{ne,y) (27) 

n odd even 

where J2odd (Seven) ^ ^^"^ restricted to odd (even) integer. Since the thermodynamical properties of the system 
are invariant under the shift q q + 2ne, one should note that g{q + 2ne,y) — g{q,y), where n is an integer. As a 
consequence, Eq. (|27|l leads to 

- =Pog(0,y)+Pe5(e,2/) (28) 
a 

where po = X^e^en (Pe ~ J2oddPn) probability that the offset charge q is an even (odd) integer multiple 

of e. The results obtained from Eq. I|27|l with the probability distribution (|26|l are displayed in Fig. El where we 
plot the phase boundary line for p^ = 1/4, 1/2, 3/4. One observes that increasing pe leads to an enlargement of the 
superconducting phase. It is worth noting that applying the SCMFT approximation with the probability distribution 
(|26|l it is possible to find exactly the same result obtained in a path-integral approach for a JJA model with diagonal 
capacitance matrices |l3j. 
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FIG. 2: Phase diagram for random offset charges with the probability distribution given by Eq. 12611 . In the plot we take 
Pe = 1/4 (dashed line), 1/2 (solid line), and 3/4 (dotted line). 

B. Random Self-Capacitances 

In this Section we shall study the finite temperature phase diagram of JJA's with uniform charge frustration q and 
random self-capacitances Cu with average Cq. Correspondingly, the average charging energy is Eq = e^/2Co. It is 
useful to define the charging energy terms Uu = SE^^ and Uq — SE^; averaging the self-consistency equation H25|l 
over the random variables Uu , the equation for the phase boundary becomes 

1 f°° P(U) 

-= / dU^g{q,U,y) (29) 
a Jo U 

where now a = 2zEj/Uo, y = SkBTc/Uo and U — Uu/Uq. Equation H29|l can be also obtained by using the path- 
integral approach Il3l|. The fimction g{q — e, U, y) is given by 

iU+y fT^lv + V 1 -(4)7/y)(«+l/2)^ 

4iy ^ ^ l-4(n+l/2)^ ^ 

g(g = e, [/, y ) , (30) 



whereas the function g{q = 0,U, y) is given by 



giq = 0, [/, y) = ^ . (31) 

1^ g-{iU/y)m^ 



If one, for instance, considers a bimodal distribution of the C/'s, then 

PiU)^p5{Ui~U) + il-p)5iU2-U), (32) 
where Ui and U2 are positive numbers. Inserting the probability distribution 132|l in Eq. I|29(l . one gets: 

- = P 77- 5(9, Ui,y) + {l~p)-^ g{q, C/2, y). (33) 

The phase boundary line given by Eq. H33() is plotted in Fig. [3] One observes that, when q = 0, the superconducting 
phase increases in comparison to the nonrandom case: this is due to the factors 1/C/i and I/U2 in Eq. (|33|l . which 
make larger the contribution of junctions with charging energies less than Uq. The increase of the superconducting 
phase is thus due to a decrease of the effective value of the charging energy. This phenomenon is largely independent 
from the specific choice of the distribution [T3| . 

It is pertinent to observe that, when q ~ e (maximum frustration induced by the external offset charges), the 
randomness does not modify considerably the phase diagram. This should be compared with the unfrustrated case 
(q = 0), where randomness sensibly affects the phase diagram. 
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(a) (b) 




1 

a 



FIG. 3: Phase diagram in the Tc — a plane for random diagonal charging energies distributed according to Eq. 13211 and uniform 
offset charge g = (a), q = e (b). Uo is the average charging energy. We take U\ = 1/2 and 172 = 3/2, while p is 1/4 (dashed 
line), 1/2 (solid line) and 3/4 (dotted line). 



IV. CONCLUDING REMARKS 

In this paper, we reviewed the use of the self-consistent mean-field theory to analyze the effects induced by offset 
charges on the finite temperature phase diagram of Josephson junction arrays. 

We reviewed, for a diagonal Coulomb interaction matrix, the explicit derivation of the equation for the phase 
boundary line between the insulating and superconducting phase. The resulting phase diagram is drawn for a uniform 
offset charge distribution q: with q = e, the superconducting phase increases with respect to g = 0, and the model 
exhibits superconductivity for all the values of a = zEj/AEc- An offset charge q = e tends therefore to decrease the 
charging energy and thus favors the superconducting behavior even for small Josephson energy Ej. 

Using a pertinent extension of the self-consistent mean-field approach, we obtained here also the phase diagram at 
finite temperature of JJA's with capacitive disorder. For a random distribution of offset charges which are integer 
multiples of e one has that the superconducting phase increases in comparison with the unfrustrated case. For 
arrays with random charging energies, the superconducting phase increases with respect to the situation in which all 
self-capacitances are equal. 

It is comforting to see that our mean-field analysis provides results which are in very good agreement with those 
obtained by Quantum Monte Carlo simulations [ij for most of the phase diagram. 
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